Uniting Bose-Einstein condensates in optical resonators 
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The relative phase of two initially independent Bose-Einstein condensates can be laser cooled 
to unite the two condensates by putting them into a ring cavity and coupling them with an 
internal Josephson junction. First, we show that this phase cooling process already appears 
within a semiclassical model. We calculate the stationary states, find regions of bistable behavior 
and suggest a Ramsey-type experiment to measure the build up of phase coherence between 
the condensates. We also study quantum effects and imperfections of the system. 
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During recent years Bose-Einstein condensates (BEC) 
of Alkali atoms and Hydrogen have been produced and 
studied extensively in the laboratory (!]]. In view of 
potential applications, such as the generation of bright 
beams of coherent matter waves (atom laser), a central 
goal has been the formation of condensates with a num- 
ber of atoms as large as possible. It is thus of particular 
interest to study a scenario where this goal is achieved by 
uniting two (or more) independently grown condensates 
to form one large single condensate. Physically speaking, 
two independently formed condensates are characterized 
by a random relative phase of their macroscopic wave 
functions. A "fusion" of two condensates thus amounts 
to locking the relative phase in a dissipative process. Be- 
low we will study such a mechanism in the context of 
optical Cavity QED @. In our schemes two condensates 
in different internal atomic states are coupled by lasers 
in a Raman configuration (internal Josephson effect |3|), 
and in addition to a lossy optical cavity, which provides 
an effective zero temperature reservoir. As a result, we 
obtain a damping mechanism for the relative condensate 
phase. A physical picture behind this cooling mechanism 
can be given by establishing a formal analogy of the dy- 
namical equations describing this "laser cooling" of the 
relative phase, and the recently discussed cavity assisted 
laser cooling of the motional degrees of freedom of atoms 
or molecules j|. 

We consider two equally large, independently produced 
BECs, where the atoms of the individual condensates are 
the same species but in two different hyperfine states 
|1) and 1 2). The elements of the one-body density ma- 
trix are given by p fei (x,x') = (tpj (x')^fc(x)) @- Here, 
^fc(x) is a bosonic field operator which annihilates a 
particle at position x and in hyperfine state |fc), and 
k,l G {1,2}. The one-body density matrix of these 
two independent condensates has vanishing off diago- 
nal elements pi2(x, x') = 0, and two dominant eigen- 
values of approximately N/2, where N is the total par- 
ticle number. The diagonal terms can then be written 
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FIG. 1: a) Experimental setup: Two initially independent 
condensates 1 and 2 are trapped in a ring cavity and coupled 
to the cavity mode c as well as to the two lasers £1 and fii. 
b) Level structure: The Raman lasers O transfer atoms from 
1 1} to 1 2) . A laser fii drives the transition between particles 
in 1 1) and an auxiliary level |3). The cavity couples the levels 
1 3} and |2) with coupling strength g c . 



as pfcfc(x,x') « Nipl(x.')(pk(n)/2, with ^fe(x) the wave 
function of the condensate in state \k). There is thus no 
phase-coherence between the two different condensates 
and any attempt to measure a relative condensate phase 
$ will produce a random result {1. Here we study the 
possibility of building up and locking the phase $ by (i) 
an internal atomic Josephson junction (J J) || [| transfer- 
ring atoms between the two initial condensates, and (ii) 
coupling the atoms to a dissipative ring cavity |Q. The 
final result is that for experimentally relevant parame- 
ters the phase $ develops towards a stationary value on 
a time scale of a few trap oscillation periods. This is re- 
flected by one dominant eigenvalue of approximately N 
of /9fej(x, x'). The corresponding condensate mode may 
be an electronic superposition state. In this sense we 
have thus joined two condensates together to form one 
larger condensate. 

The internal JJ || is realized by coupling the two hy- 
perfine states {|1), 1 2)} by a Raman transition with de- 
tuning A and effective Rabi frequency Q (cf. Fig. [l]). 
We assume f2 to be real, positive, and x independent, 
i.e. there is negligible momentum transfer to the conden- 
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sates due to the Raman lasers. In second quantized form 
the Hamiltonian Pbec = Pi + H 2 + Pl H is given by 
(with h = 1) 
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where ipifi = ^1,2 (x), V = V(x) the trapping poten- 
tial; we have suppressed x for notational convenience. 
We assume the trapping potential to be harmonic with 
frequency u> and the mass of the atoms to be m. Two- 
particle interactions are characterized by Ukk — ^a^jm 
where a k is the s-wave scattering length of atoms in con- 
densate k and by U12 = U21 = ina^/m with a\% the 
interspecies s-wave scattering length. We restrict our- 
selves to the case where a\a-2 > af 2 , i-e. to two stable 
strongly overlapping condensates || . 

The ring cavity mode couples the state |2) to an aux- 
iliary excited atomic state |3), with detuning S c and cou- 
pling strength <? c (x), which decays at a rate re. We also 
add a classical laser field driving the |1) to |3) transition 
with Rabi frequency fii(x) and detuning Si (cf. Fig. [I]). 
Adiabatically eliminating the internal state |3) (which re- 
quires Si » Oi(x) and S c » g c (x)\C\, with |C| the square 
root of the number of photons in the cavity) we obtain 
the following master equation with H = Pbec + He'- 
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where c is the annihilation operator of a photon in the 
cavity mode, v = S c — Si — A is the effective cavity de- 
tuning, and 



He = dx g(x)cip\ip2 + h.c. + 
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The effective coupling of the condensates to the cavity 
is given by g(x) = Qi(x)g c (x)/25i. In a ring cavity the 
spatial dependences of fii(x) and <? c (x) can be made to 
almost cancel each other. There is thus negligible mo- 
mentum transfer to the atoms due to the coupling to the 
cavity, and we may set <?(x) = g, where g is x indepen- 
dent, real, and positive. 

We derive equations of motion for the operators tpk (x) 
and c. First we study laser cooling of the phase $ in a 
semiclassical model before returning to the full quantum 
model below. In the semiclassical model we assume N 3> 
1 and replace the operators ?pk (x) and c by the c- numbers 
*fe(x) and C, respectively. We define C\ = C, C 2 = C*, 
denote the Kronecker delta by Ski and find 
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In terms of this semiclassical framework the initial inde- 
pendence of the two condensates is modelled by a ran- 
dom initial phase, i.e. pki(x,x') » (4 , *(x')Wfc(x)) c i. Here 
(. ..) c ] is the stochastic average over relative phases $ 
which are initially equally distributed between [0,27r]. 
Apart from the relative phase the initial conditions for 
solving Eqs. (Q) are assumed to be identical. 

For analytical calculations our numerical studies 
(see below) suggest the two mode ansatz (with <p(x) 
a fixed and normalized wave function) ^(x) = 
ViVfc</5(x)exp(i$fe), where Nk is the expectation value of 
the number of particles in state \k), and is an x inde- 
pendent total phase. Within this semiclassical two mode 
model we find: 
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where we have defined $ = $ x - $ 2 , P$ = (N 2 - Ni)/2 
and a = (gC + £1/2) exp(— i$). The chemical potentials 
p,k are defined by 
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To make contact with the resistively shunted junction 
model ||, often used in describing JJs we set un = u 22 , 
and assume N > 2|P$| and NU > {CI, gC}, which 
can be enforced by adjusting the laser parameters cor- 
respondingly. Here U = («u — U12) f dx|</j(x)| 4 . We 
further simplify Eqs. (j^) by adiabatically eliminating 
the cavity mode C to second order, valid when re 3> 
{gyN,UN}, and derive an intuitively appealing damp- 
ing equation for a fictitious particle moving along a co- 
ordinate $ (with M = 1/2U), 
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The particle moves in a potential V($) = CIN cos($)/2 + 
Prf$ where Fd = g 2 N 2 n/{n 2 + 4is 2 ) is a constant 
force. This potential has minima at $ = (2n + 1)7T — 
arcsin(r2 c /ri) with integer n for fl > Sl c = 2Fg/N. 
Close to a minimum V(&) is harmonic with frequency 
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The motion of the particle 



in V($) is damped with a friction coefficient a r = 
vg 2 N 2 n/2 (k 2 /4 + v 2 ) 2 . From Eq. (§) we easily find the 
damping time scale of the relative phase which is given 
by t = 2M/a r . Note that cooling of the phase towards 
a minimum of the potential occurs only for a r > which 
requires a cavity detuning of v > and that both Fd and 
a r are induced by the cavity damping rate re. 
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If NU < {Q,gC} Eq. (0) no longer holds. However, 
from Eqs. @ we can still find an approximate expression 
for r which has a minimum at Q = Q m = [{Av 2 — n 2 + 
2V« 4 + Av 2 n 2 + 16i/ 4 )/12] 1 / 2 . This case corresponds to 
the Rabi-oscillation limit ||. 

We now determine the stationary solutions of our semi- 
classical two mode model [Eqs. (^)] which are charac- 
terized by C = 0, P$ = and by $ = 0. Writing 
C = \C\ exp(i$ c ) we find the following conditions: 



2. 9v /A 1 A 2 sin($-$ c ) =k\C\, 
- 5 V^V2Cos($-$ c ) =v\C\, 
-Ov/ivTiVasin^ =k\C\ 2 , 

(Nx - iV 2 )Re{a}. 
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We set the macroscopic wave function <ys(x) = <ySo( x ) 
where v?o( x ) is found from numerically solving Eqs. ( pla| ) 
for Ni = N2 with all the lasers turned off (O = g = 
A = 0) for the ground state. Thus we find ^2 — ^x = 
U(N2 — Nx). As shown in Fig. || we obtain one sta- 
ble (unstable) branch A (B) of stationary solutions with 
Ni = N 2 for fl > fl c . For fl < fl T: where fl 2 , = 
il 2 c + [NU + Nvg 2 /(v 2 + k 2 /A)] 2 , there exist stationary 
solutions with Nx ^ N2 where the solutions on branch C 
(D) with Ni < N/2 (Nx > N/2) are stable (unstable). 
The solutions of branches A and B can also be found from 
the Josephson model Eq. (Q) while branches C and D are 
beyond the range of validity of Eq. (0). The stability of 
the stationary solutions is checked by a linear stability 
analysis of Eqs. (||). We also numerically evolve Eqs. (^) 
with initial using the two mode ansatz where Nx, N2, $, 
are found from Eqs. (||) and <p(x.) = <po(x-)- In numeri- 
cally solving Eq. (^) we restrict ourselves to one spatial 
dimension ||. For solutions on branch A these initial 
conditions are exact stationary solutions of Eqs. (|j) and 
we find them to be numerically stable. The phase $ 
always evolves towards its stationary value. Therefore 
solutions on branch A are best suited for uniting two 
condensates. On branch C with Nx 7^ N2 which is less 
interesting for the purpose of laser cooling of the phase 
5> the validity of the two mode ansatz depends on the 
strength of the two particle interaction. For UN < u> 
we find good agreement between the numerics and our 
analytical results while for UN > lo instabilities in the 
evolution of Eqs. ([|) emerge which do not appear in the 
linear stability analysis of Eqs. (||). 

A possible experimental scenario to monitor the lock- 
ing of the relative phase is a Ramsey experiment. Af- 
ter time t cooling is stopped and a tt/2 pulse with phase 
§r is applied to the condensates ||. This pulse effec- 
tively transforms the relative phase $ into a difference 
in occupation numbers. For a well defined phase by 
varying we expect to observe maximum fringe visi- 
bility v = 1, whereas for an undefined, effectively ran- 
dom phase v = 0. In terms of the one particle den- 
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FIG. 2: Mean particle number iVi in state |1) [a)] and phase 
$ [b)] of the condensate in the stationary state against Rabi- 
frequency fl. The solid (dashed) curves correspond to stable 
(unstable) stationary solutions. The dotted vertical lines give 
the critical Rabi frequencies Q c and Qt as defined in the text. 
The parameters are UN = 0.213a;, A = 0, \fNg = 3.87a;, 
k = 40a> and v = 20a;, corresponding to condensates of total 
N = 6000 Rb atoms with m = 1.44 • 10~ 25 kg and scattering 
lengths of ai = 02 = 5.45nm, a\2 = 5.24nm trapped in a 
harmonic potential with frequency u> = 200Hz/27r. 



sity matrix pfci(x, x') the visibility is given by v(x) = 
2|pi 2 (x,x)|/[pn(x,x) +p 2 2(x,x)] at position x. Fig- 
ure ||a shows the visibility v = u(x = 0) in the center 
of a one dimensional trap found by numerically solving 
the semiclassical model [Eqs. (^)] against the cooling time 
t. As expected v does not go to 1 in the bistable region 
while it tends to 1 for > f2j\ The results agree very 
well with the semiclassical two mode model Eqs. (||). 

So far we have studied a semiclassical version of a 
Josephson-cavity model, where the atoms and the cavity 
mode are described by c-number fields. We will investi- 
gate now quantum effects within a two-mode model of a 
J J ||. This is obtained from the Hamiltonian (|l|) with 
the replacement ^fc(x) = ^(^(x) where bk are destruc- 
tion operators of a particle in electronic state |fc), and 
</?(x) a (fixed) spatial mode function. We define an op- 
erator describing the imbalance in the particle numbers 
P0 = (b\b 2 — b\bx)/2, and its canonical conjugate —<j>, 
where (f> has the meaning of a relative phase operator Q . 
Projection on the eigenstates of the phase operator, \cf>), 
gives the quantum phase distribution p(tfi) = tr{|0)((/)|p}. 
We note that the eigenstates of the phase operator \cf>) ~ 
Em e- im *| N/2 - to) 1 1 N/2 + m) 2 are entangled states of 
particles in the first and second internal state jjj. For 
small particle numbers and low occupation of the cav- 
ity mode, the master equation (||) can be solved directly 
by quantum Monte Carlo techniques, and the time evo- 
lution of the phase distribution can thus be calculated 
exactly. As an example, Fig. ||b shows the evolution of 
the phase distribution p(4>) for N — 50 atoms, Uxx = u 22, 
and initially independent condensates, i.e. p(4>) = 1/2tt, 
obtained from a simulation. During the cooling process 
p{4>) builds up around the semiclassical stationary solu- 
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tion for <E> on the time scale predicted by the semiclas- 
sical model, and one finds excellent agreement between 
the Monte Carlo simulation and the semiclassical model, 
even for a small number of particles. 

The quantum limits of the steady state distribution 
of the phase can be discussed within a quantum phase 
model. For UN ^> {ft, gc} the Hamiltonian simplifies to 
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which describes the motion of a fictitious quantum par- 
ticle in a cose/) phase-potential coupled to a (damped) 
harmonic oscillator representing the cavity mode. In this 
picture, cavity assisted phase locking can be understood 
as "cooling of the motion" of this fictitious particle in the 
phase potential due to coupling to a damped harmonic 
oscillator. In fact, the master equation (||) with Hamil- 
tonian (|^) is identical to the master equation which have 
been derived for laser cooling of single atoms or molecules 
in a high-Q cavity moving in a trapping potential [Q. 
Analytical expressions for the width of the steady state 
phase distributions can be given by expanding the Hamil- 
tonian (^) around the semiclassical steady state solutions 
to second order, and solving for the variance of <f> and p^ 
in the stationary state. We find (^] 

2 _ k 4 + Ak 2 (uj 2 + 2v 2 )~16fvuj+16v 2 (u; 2 + v 2 ) 
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where g = Ngy/U/u. In the limit of cavity width much 
larger than the oscillation frequency in the phase po- 
tential, k 3> {Co,g}, the minimum uncertainty of the 
phase variance is obtained for the detuning v = k/2, 
and is given by A(f> 2 Uk/2lo 2 ^> U/ui. On the other 
hand, for small cavity damping, {n,g} -C {v, Co}, we find 
Acj) 2 w U/ui and Ap^ « Co/AU for v = ui. By analogy 
with laser cooling we can identify the first limit with the 
Doppler regime, where a semiclassical description is valid 
and the final temperature is kT — k/2, whereas the sec- 
ond case corresponds to the sideband cooling limit, where 
the particle is cooled to the vibrational ground state, i.e. 
we have a minimum uncertainty state (T = 0). 

The above results were obtained for space independent 
couplings g(x) = g and fi(x) = fl. However, due to im- 
perfections absorption and emission of photons may lead 
to a small momentum transfer to the condensates result- 
ing in spatial perturbations to the positional condensate 
mode <£>(x). A detailed analysis shows that these per- 
turbations are damped away with a cooling mechanism 
similar to the phase cooling described above . 

In summary, we have shown that by coupling two con- 
densates to a dissipative ring cavity we can "laser cool" 




FIG. 3: a) Visibility v in the Ramsey experiment against cool- 
ing time t for the parameters as in Fig. Be Q, = 15lu (solid line), 
fl = 0.55oj (dashed) where two stable stationary solutions ex- 
ist, b) Time evolution of the phase distribution p(<fi) from a 
Monte-Carlo simulation. The dashed vertical lines indicate 
the semiclassical solution in the stable stationary state. For 
parameters see Fig. |^, f2 = 40a; and TV = 50 particles. 



the relative phase between two initially independent con- 
densates without loss of particles. 
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